Thermal Relic Abundances of Particles with Velocity-Dependent Interactions 
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We reexamine the evolution of thermal relic particle abundances for the case where the interaction 
rate depends on the particle velocities. For the case of Sommerfeld enhancement, we show that the 
standard analytic approximation, modified in a straightforward way, provides an estimate of the 
relic particle abundance that is accurate to within 10% (in comparison to < 1 % error for the non- 
Sommerfeld-enhanced case). We examine the effect of kinetic decoupling on relic particle abundances 
when the interaction rate depends on the velocity. For the case of pure p-wave annihilation, the 
effect of kinetic decoupling is an increase in the relic abundance, but the effect is negligible when 
the kinetic decoupling temperature is much less than the chemical decoupling temperature. For the 
case of Sommerfeld-enhanced s-wave annihilations, after kinetic decoupling occurs, annihilations 
continue to change the particle abundance down to arbitrarily low temperatures, until either matter 
domination begins or the Sommerfeld effect cuts off. We derive analytic approximations to give the 
final relic particle abundances for both of these cases. 
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I. INTRODUCTION 



The calculation of the thermal evolution of particle 
abundances in the early universe represents one of the 
earliest and most fundamental applications of particle 
physics to cosmology [H4l2l]. The most important appli- 
cation of this calculation is the determination of the relic 
dark matter abundance. For the simplest case, that of 
s-wave annihilation, one finds that the final abundance 
is essentially independent of the mass, and the observed 
dark matter abundance can be achieved with an annihi- 
lation rate of [l3[ approximately (cry) ~ 3 x 10~ 26 cm 3 
sec -1 , where (av) is the thermally averaged product of 
the cross section and relative velocity of the annihilating 
dark matter particles. This result is quite interesting, as 
it suggests that physics at the electroweak scale may be 
responsible for the dark matter. 

Recently an interesting twist has emerged in this calcu- 
lation. Motivated by a desire to explain various anoma- 
lous astrophysical backgrounds, a number of investigators 
have examined the possibility that dark matter annihi- 
lation involves a Sommerfeld enhancement, which pro- 
vides an additional factor of 1/v in the dark matter an- 
nihilation cross section [l4| [l5[ . The effect of this Som- 
merfeld enhancement on the thermal relic abundances 
has been discussed for specific models in [l6l - [l9j . and 
treated more generally by Kamionkowski and Profumo 
[l4| and Arkani-Hamed et al. [Hj]. Note that the latter 
two papers reached opposite conclusions regarding the 
effect of Sommerfeld-enhanced annihilations on the relic 
abundance. Kamionkowsk and Profumo estimated a sig- 



nificant suppression, while Arkani-Hamed et al. argued 
for a very small effect. However, these two conclusions 
are not actually inconsistent, because they rely on differ- 
ent assumptions regarding the strength of the coupling 
that induces the Sommerfeld enhancement. We will ex- 
amine both limiting cases in our discussion below. 
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In this paper, we consider several new aspects of veloc- 
ity dependent interactions, including both Sommerfeld- 
enhanced s-wave annihilations, and pure p-wave annihi- 
lations. In the next section, we first consider the an- 
alytic approximation of 0, modified for the case 
of Sommerfeld-enhanced s-wave annihilations, and show 
that it provides a good approximation to the relic abun- 
dances in this case. We then consider the effects of ki- 
netic decoupling, which increases the rate at which the 
relic particle temperature declines, and thereby modi- 
fies the abundance when the annilation rate has a veloc- 
ity dependence. We provide estimates of this effect for 
both pure p-wave annihilation and Sommerfeld-enhanced 
s-wave annihilation. The latter is a significantly larger 
effect; we find that in this case, annihilations continue 
to decrease the relic abundance down to arbitrarily late 
times, until the abundance freezes out either at the on- 
set of matter domination, or when the Sommerfeld effect 
itself cuts off. Although calculations of this sort can al- 
ways be done numerically for any particular model of a 
relic particle, it is useful to derive such analytic estimates, 
since they can be applied to arbitrary models, and can 
provide qualitative insight into the behavior of such mod- 
els. While interest in the Sommerfeld enhancement has 
been spurred by recent astrophysical observations, our 
discussion here is intended to be as general as possible. 
Our results are discussed in Sec. III. 



2 



II. CALCULATION OF RELIC ABUNDANCES 

A. Sommerfeld-enhanced s-wave annihilation 

Recall first the standard formalism for thermal parti- 
cle abundances in the early universe [9T4TTj]. Let n be 
the number density of a relic particle and n eq be its 
thermal equilibrium number density. Then 



^ + 3Hn = -{av)(n 2 -nl q ), 



(1) 



where H is the Hubble parameter. To eliminate the ex- 
pansion term, we express the number density in terms of 
Y = n/s, where s is the total entropy of the universe, and 
we change the indendent variable to x = m/T. Further, 
following [H, [l(| , we parametrize the cross-section as 



(av) = (To x 



(2) 



where n = corresponds to s-wave annilation, n = 1 for 
p-w&ve annihilation, and so on. Note that [ll[ provides a 
more sophisticated treatment of (crv), but at the level of 
accuracy we are interested in here, equation ([2]) will be 
sufficient. For all of the specific cases examined here, we 
take Co = 3 x 10 -26 cm 3 s _1 . Then equation ([IJ becomes 



dY 
dx 



Y 2 ), 



-Xx- n - 2 (Y 2 
where the constant A is given by 

A = yj tt /45(5»s / 'gl /2 )m Pl m x cr , 



(3) 



(4) 



with mpi = 1/y/G. Here g* is the effective number of 
relativistic degrees of freedom in the universe, defined 
by the requirement that the energy density in relativistic 
particles is pn = (tt 2 /30)g*T 4 , while g* s is defined in 
terms of the entropy density s as s = (2ir 2 / 45) <?* S T 3 . For 
the cases we examine here, it is accurate to take <?* s g», 
and is given by 

.9, = 106.75 T > 175 GeV (5) 

.9, = 96.25 175 GeV > T > 80 GeV (6) 

.9, = 86.25 80 GeV > T > 4 GeV (7) 

g* = 75.75 4 GeV > T > 150 MeV (8) 

g* = 17.25 150 MeV > T > 20 MeV (9) 

.9, = 10.75 T < 20 MeV (10) 

For the case of interest here, the relic particles can 
assumed to be nonrelativistic, so that Y eq is well- 
approximated by Maxwell-Boltzmann statistics: 



Y EQ = .U5(g x /g*)x^e 



3/2^* = nrr 3/2 t 



(11) 



where g x is the number of degrees of freedom of the x 
particle. 



At early times, the relic particle is in thermal equilib- 
rium, so that its abundance tracks the equilibrium abun- 
dance, but at late times the abundance freezes out to a 
constant value. This argument can be made more explicit 
by defining the quantity A = Y 
which is given by 



Y eq , the evolution of 



dA 

dx 



dY eq 
dx 



Xx- n - 2 A(2Y eq + A) 



(12) 



The approximation in [9( and [10( amounts to setting the 
right-hand side of equation (|12[) to zero up to xj, the 
value of x at which the abundance freezes out, and then 
integrating equation (|12[) for x > Xf with the assump- 
tion that both Y eq and dY eq /dx are negligible. One then 
obtains @, OH 

x f = ln[(n + l)a\] - (n + l/2)ln[ln[(n + l)oA]], (13) 



and the final value of Y is 

3.79(n+ l)x n f 



n+l 



Y„ 



(g* s I gl /2 )mpim x (TQ 



(14) 



This approximation agrees with the exact integration of 
the Boltzmann equation to within a few percent. For 
the case of s-wave annihilations, the evolution of A is 
compared to the approximate evolution in Fig. 1 for a 
500 GeV particle. 

Now consider what happens for s-wave annihilations 
that are Sommerfeld enhanced. Sommerfeld enhance- 
ment arises from a long-range attractive force due to a 
light force carrier (f>. In the limit where — > 0, the 
annihilation cross-section is enhanced by the factor [l5| 



S = 



ixa/v 



\ g—na/v ' 



(15) 



where v is the velocity of the annihilating particles, and 
Aira is the square of the coupling. 

Clearly, the magnitude of the Sommerfeld enhance- 
ment depends on the value of a. We illustrate this ef- 
fect in Fig. 1, showing how the evolution of the particle 
abundance depends on a. Clearly, for a ;$ 0.01, the ef- 
fect of the Sommerfeld enhancement on the final relic 
particle abundance is negligible. On the other hand, for 
ita/v 3> 1, equation (1151) reduces to a l/v enhancement 
in the annihilation cross-section. Fig. 1 shows that this 
limit is achieved for a 



15] 



as- 



0.3. Note that [14| and 
sumed these opposite limiting behaviors in deriving their 
estimates of the effect on the thermal relic abundance. 
The case 0.01 ^, a ^ 0.3, gives an intermediate regime 
displayed in Fig. 1. Since l/v enhancement provides one 
set of limiting behaviors, we will assume a simple l/v 
enhancement in what follows. 

Note that the effect is more complex if one does not 
assume — > 0; in this case, the production of bound 
states results in resonant enhancement of the annihilation 
rate, while at the same time the Sommerfeld enhance- 
ment cuts off for v < 771^ /m x 15j . We will consider only 
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no Sommerfeld enhancement 




Sommerfeld a -.01 



Sommerfeld a = I 



1 Soninicrlcld a =.3 



FIG. 1: The evolution of relic particle densities for the case 
of Sommerfeld-enhanced s-wave annihilations of a 500 GeV 
mass particle with ao — 3 x 10 -26 cm 3 s^ 1 as a function of 
the indicated value of the coupling a. Top and bottom curves 
correspond to the limiting cases of no Sommerfeld enhance- 
ment, and 1/v enhancement. Horizontal lines are the analytic 
estimates for the final relic abundances in these two cases (i.e. 
n = and n = —1/2, respectively, in equations 1 131 and 1 14|) . 
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FIG. 2: A comparison between the analytic approximation for 
the ratio in relic abundance for s-wave annihilation with and 
without a 1/v Sommerfeld enhancement to the correspond- 
ing numerical results over several orders of magnitude of the 
parameter A. 



the case — > 0, but will discuss these other effects 
qualitatively later. 

A nonrelativistic particle has a velocity that scales as 
(v 2 ) oc T x . As long as the particle is in thermal equi- 
librium, T x = T. Hence the effect of Sommerfeld an- 
nihilation is to modify equation ([3]) (for n — 0) to the 
form 

d J- = -\ x -*l\Y^Y 2 q ). (16) 

It would appear, then, that the relic abundance in this 
case is well-approximated by the standard freeze-out 
abundance for the case n = —1/2 in Eqs. JT3J) and (ITU) . 
Indeed, this was the assumption made in Ref. [TiJ . How- 
ever, it is not a priori obvious that equations (|13[) and 
(1141) can be accurately applied in the regime n < 0, since 
they have been numerically tested only in the regime 
n > 0, and the freeze-out process becomes progressively 
less "sharp" as n decreases. In Fig. 1, we integrate the 
Boltzmann equation for the same set of parameters, but 
with a 1/v enhancement in the annihilation rate. The 
final abundance is reasonably well-approximated by the 
n = —1/2 analytic approximation, but the agreement 
with the exact numerical results is not quite as good as 
for the s-wave case without the Sommerfeld effect. 

This result allows us to estimate the ratio between the 
abundance in the presence of Sommerfeld enhancement, 
s-wave abundance without Som- 
merfeld enhancement, Y^. We obtain 

■^oo 2 X j 

1 v /ln(aA/2) 

21n(aA) - (1/2) lnln(aA) ' 1 ' 



Here xjsom is the value of Xf when Sommerfeld enhance- 
ment is included. Taking x fgp u ~ a;/o yields the abun- 
dance estimate given in [l4|. Our results confirm that 
the change in the Xf is indeed very small. For the stan- 
dard dark matter freeze-out value of Xf = 20, we find 
that XfsoM is larger by only a few percent, while equa- 
tion CE7D gives Y ( ^ OM /Y OD - 1/10. Both of these results 
are confirmed numerically in Fig. 1. In Fig. 2, we com- 
pare the estimate given by equation dTTl) to numerical 
results. The analytic estimate differs from the numerical 
result by about 10%. In contrast, the analytic abundance 
estimate for s-wave annihilation without Sommerfeld en- 
hancement is accurate to within 1% M. 



B. Effects of Kinetic Decoupling 

The results in the previous section assume that the 
temperature of the relic annihilating particle tracks the 
background radiation temperature. This will be true as 
long as the annihilating particle remains in thermal equi- 
librium with the radiation background. However, once 
the particle drops out of thermal equilibrium, we expect 
its temperature to scale as T x oc 1/R 2 , where R is the 
scale factor, while the radiation temperature scales as 
T oc 1/R. 

Thus, we need to make the standard distinction be- 
tween chemical equilibrium and kinetic equilibrium. The 
freeze-out process we have discussed above is actually the 
process by which the annihilating relic particle drops out 
of chemical equilibrium, so that number-changing inter- 
actions are no longer effective, and the particle's comov- 
ing number density becomes constant. However, even 
after dropping out of chemical equilibrium, the relic par- 
ticle will, in general remain in kinetic equilibrium as it 
continues to scatter off of relativistic standard model par- 
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tides which are in local thermal equilibrium with the ra- 
diation background. As long as the particle is in kinetic 
equilibrium, its temperature tracks the background ra- 
diation temperature. Finally, at some kinetic decoupling 
temperature, T^, the scattering interactions are no longer 
sufficient to maintain kinetic equilibrium, and the tem- 
perature of the particle decreases as 1/R 2 rather than 
1/R. (For a recent discussion, see, e.g., [20I4231 ] V 

The precise temperature at which kinetic decoupling 
occurs is dependent on the model for the relic particle of 
interest. For instance, in the set of supersymmetric mod- 
els examined in [23], T k /T f ~ f (T 1 - 1CT 3 . Since we wish 
to keep our discussion as general as possible, we will take 
Tk/Tf as a free parameter, subject only to the constraint 
that Tk < Tj, since number-changing interactions also 
maintain kinetic equilibrium. We also make the approxi- 
mation that the particle drops out of kinetic equilibrium 
instantaneously at Tk ; this is a reasonable approximation 
[23j . With these assumptions, the relation between 



T x and T is 



T x = T 2 /T k 



(18) 



The change in the evolution of T x brought about by 
kinetic decoupling changes the velocity evolution of the 

1/2 

annihilating particles, since v cx T£ for nonrclativistic 
particles. Since the standard s-wave annihilation cross 
section, (av), is independent of T x (or equivalently, v x ), 
kinetic decoupling has no effect in this case. The same is 
not true for p-wave annihilation, for which (av) cx T x , or 
for Sommerfeld-enhanced s-wave annihilation, for which 
(av) cx T x 1 . The reverse reactions (which create x) can 
be neglected during the era following kinetic decoupling, 
since Tk < Tf. Thus, the Boltzmann equation following 
kinetic decoupling for p-wave annihilation becomes 



%- = -Xx k x 
ax 



-4y2 



(19) 



where we define the constant Xk = ni/Tk- For 
Sommerfeld-enhanced s-wave annihilation, we obtain: 



~ ~ XXk X Y ■ 



(20) 



The effect of kinetic decoupling on the final relic abun- 
dances is easy to estimate. Recall that equation ([14)) 
is derived by integrating the annihilation portion of the 
Boltzmann equation (only) from x = Xf to 00 0, ITol . 
Replacing this integration by an integration from Xf to 
Xk , and then integrating equations (fT9)) and ([20)) from Xk 
to 00 should provide the correct estimate of the change 
in the final relic abundance. 

For p-wave annihilation, we obtain the ratio between 
the final abundance in the presence of kinetic decoupling, 
Yk?\ and the abundance in the limit where the particle 
stays in kinetic equilibrium to an arbitrarily low temper- 
ature, Yoo. This ratio is 



Y- 



(*) 



1 



(l/3)(T k /T f y 



(21) 




FIG. 3: The effect of kinetic decoupling on the evolution of 
the relic particle abundance for the case of p-wave annihilation 
of a 500 GeV mass particle with <tq = 3 x 10~ 26 cm 3 s _1 . 
Horizontal line gives our analytic estimate of the final relic 
abundance. 



We see that the effect of kinetic decoupling is to in- 
crease the final relic abundance for the case of p-wave 
annihilations. This easy to understand, since the anni- 
hilation rate in this case scales as T x , so a more rapid 
decrease in T x due to kinetic decoupling leads to fewer 
relic annihilations after freeze-out, and so a larger relic 
abundance. The effect, however, rapidly becomes irrel- 
evant for Tk/T f « 1. For example, for T k /T f = 1/2, 
the result is a 9% increase in the relic abundance. For 
Tk/Tf < 0.1, the increase in the relic abundance is less 
than 0.3%. A numerical calculation of this effect is illus- 
trated in Fig. 3. Since the effect of kinetic decoupling 
becomes significant only for values of Tk /Tf that are im- 
plausibly large, it is unlikely to be of much importance 
for p-wave annihilation. 

The effect of kinetic decoupling is much more strik- 
ing for the case of Sommerfeld-enhanced s-wave anni- 
hilations. In this case, an integration of equation ([20]) 
shows that annihilations never terminate: after kinetic 
decoupling, Y ~ 1/ In x. (This case has previously been 
discussed briefly in [9|, LL5|). However, this process will 
eventually be cut off by one of two possibilities. First, as 
noted earlier, Sommerfeld enhancement saturates once 
the velocity drops to v ~ m c f > /m x , at which point nor- 
mal s-wave annihilations resume. Second, our calculation 
holds only for the radiation-dominated case, and freeze- 
out will occur rapidly once matter domination begins. 
Let T cuto ff be the radiation temperature at which the 
Sommerfeld effect cuts off or matter domination begins, 
whichever is larger. Then we can again integrate the 
equations governing particle annihilation from Xf to Xk 
with T x = T, and from Xk to x cu toff with T x = T 2 /Tk, 
where freeze-out then occurs with negligible further an- 
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no Sommerfeld enhancement 



Sommerfeld enhancement 

No kinetic decoupling 

— i- ---- 

r, 'i'v" 



10 10" 10 10 10 10 

x = m/T 

FIG. 4: The effect of kinetic decoupling on the evolution of 
the relic particle abundance for the case of s-wave annihilation 
for a 500 GeV mass particle with oo = 3 x 10~ 26 cm 3 s~\ in 
the limit where the Sommerfeld enhancement scales as 1/v. 
Horizontal lines give our analytic estimates of the final relic 
abundances. 



When kinetic decoupling occurs, it affects the 
relic abundances for both p-wave annihilations and 



no Sommerfeld enhancement no kinetic decoupling 




Tf = T t \ 



Sommerfeld enhancement a = .01 



10 10 

x = m/T 



FIG. 5: As Fig. 4, for Sommerfeld-enhancement coupling of 
a = 0.01, a value for which the Sommerfeld effect by itself is 
negligible without kinetic decoupling. Note the strong effect 
of kinetic decoupling upon the relic particle abundances. 



nihilations at T cuto ff- We find 

^- = { Tf /T k f/\^f k -l + l -\n{T k /T cutofS ))-\ 

The effect of kinetic decoupling with Sommerfeld- 
enhanced annihilations is illustrated numerically in Figs. 
4 and 5. In Fig. 4, we show the evolution of the parti- 
cle abundance for the case we have just considered (1/v 
enhancement), while Fig. 5 shows the case a = 0.01 (of 
course, our analytic estimate, equation (|22[) , does not ap- 
ply in the latter case.) Fig. 5 illustrates the fact that a 
value of the coupling for Sommerfeld enhancement can 
be small enough to produce a negligible change in the 
relic abundance without kinetic decoupling, but it can 
have a large effect once kinetic decoupling occurs. 

III. DISCUSSION 

We have confirmed that the standard analytic approx- 
imation for the relic particle abundances can be applied, 
with the appropriate modification, to the case of s-wave 
relic abundances in the presence of a Sommerfeld en- 
hanced interaction, although the error in applying this 
approximation to the case of Sommerfeld-enhanced s- 
wave annihilations (~ 10%) is significantly larger than 
in the s-wave case without Sommerfeld enhancement (< 
1%). We have also determined the range of the coupling 
a over which Sommerfeld annihilation can be either ne- 
glected in the calculation of relic densities (as suggested 
in [Hj]) or treated purely as a 1/v enhancement to the 
annihilation rate (as in |l4|). 



Sommerfeld-enhanced s-wave annihilations. In the for- 
mer case, the effect is generally very small unless kinetic 
decoupling occurs at nearly the same epoch as chemi- 
cal decoupling. For Sommerfeld-enhanced s-wave decou- 
pling, the effect is quite large, and we have provided an 
analytic estimate of this effect. 

Finally, we note that another, quite different mecha- 
nism to produce a velocity-dependent cross section is for 
a pole to lie near twice the mass of the annihilating par- 
ticle [l2j]. The effect is most striking when the pole lies 
slightly below twice the particle mass [Hj]- In this case, 
just as for Sommerfeld-enhanced annihilation following 
kinetic decoupling, the annihilations do not freeze out 
until the velocity drops below a cut-off scale in the model. 
Since the relic abundance in his model is set by this cut- 
off scale, one would not expect a large change in the final 
relic abundance if the annihilating particles also kinet- 
ically decoupled. However, a more detailed calculation 
likes outside the scope of this paper. 
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